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Abstract 

We study the asymptotic behavior of the Hopf characteristic function of fractals 
and chaotic dynamical systems in the limit of large argument. The small argument 
behavior is determined by the moments, since the characteristic function is defined 
as their generating function. Less well known is that the large argument behavior 



* zach@het .brown.edu 
tcengizpehlevan@fas.harvard.edu 

* gerry@het .brown.edu 



is related to the fractal dimension. While this relation has been discussed in the 
literature, there has been very little in the way of explicit calculation. We attempt 
to fill this gap, with explicit calculations for the generalized Cantor set and the 
Lorenz attractor. In the case of the generalized Cantor set, we define a parameter 
characterizing the asymptotics which we show corresponds exactly to the known 
fractal dimension. The Hopf characteristic function of the Lorenz attractor is 
computed numerically, obtaining results which are consistent with Hausdorff or 
correlation dimension, albeit too crude to distinguish between them. 
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Because of the impossibility of predicting the detailed behavior of 
chaotic dynamical systems over long times, one is generally only inter- 
ested in their geometric and statistical properties. Nevertheless, the 
predominant means of computing these properties is by direct numer- 
ical simulation of a trajectory which is subsequently analyzed like ex- 
perimental data. While this is a viable approach, it is not formulated 
in terms of the quantities of interests and therefore yields little insight. 
Moreover, direct simulation is probably not the most efficient approach, 
since the simulated trajectory contains much more information than can 
be reliably calculated. Fortunately, the equations of motion yield equa- 
tions for statistical quantities, such as the Hopf characteristic function 
Z(J), which might be soluble without the need for direct simulation of 
a chaotic trajectory. 

Associated with a chaotic trajectory is an invariant measure // which 
determines the frequency with which the system may be found in a 
region of phase space. The Hopf characteristic function is the Fourier 
Stieltjes transform of this measure, Z(J) = f dp exp(zJ ■ x) The small J 
behavior of the Hopf characteristic function determines the equal time 
correlation functions of a chaotic attractor, via the Taylor Maclaurin 
series expansion in J. 

In this article, we will study the large J asymptotics of the char- 
acteristic function. While small J behavior is related to the correla- 
tion functions, the large J behavior is related to the geometry of the 
attractor, specifically its dimensionality. The geometric structure of 
chaotic attractors, and their fractal cousins, is bizarre to the unini- 
tiated. Among the strange but commonly observed features is non- 
integral dimensionality, and fine structure at all scales. With these 
features comes the attendant impossibility of defining a finite distribu- 
tion function on the attractor or fractal. The absence of a distribution 
function means that the measure dp can not be written as p(x)dx for 
any function p(x). In this case the Fourier transform of the charac- 
teristic function, J dJZ(J)exp(iJ ■ x), can not converge. The associated 
constraint on the asymptotic behavior of Z underlies quantitative state- 
ments which can be made relating the dimensionality to the asymptotics 
of the characteristic function. 

The relationship between fractal measures and the asymptotics of 
their Fourier-Stieltjes transform has been discussed in the literature 
[9]. however there has been very little in the way of explicit calculation. 
The intent of this article is to fill this gap. The Hopf characteristic 
function can be defined for a fractal, with asymptotic properties which 
are often readily computable due to the self similarity which is a hall- 
mark of fractals. We will determine asymptotics properties of Z(J) for 
the generalized Cantor set, showing the correspondence with the known 
fractal dimension. Our arguments can be readily generalized to other 
fractals. 

Determing the asymptotics of the Hopf function of a chaotic dynam- 
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ical system is a much more difficult, and mostly unsolved, problem. We 
also attempt to determine the asymptotics of the characteristic func- 
tion of the Lorenz attractor numerically. The results for the Lorenz 
attractor are crude, being limited to large but finite J, yet seem to 
be consistent with the known dimensionality. It is our hope that this 
crude approach may be supplanted by one in which the asymptotics is 
determined from the differential equations which the Hopf characteris- 
tic function satisfies, subject to the appropriate boundary conditions. 

1 Introduction 

Chaotic systems are generally studied by direct numerical simulation, with statis- 
tical (e.g. moments) and geometric (e.g. dimension) information determined from 
a time series. However, because of senstitivity to initial conditions, the detailed 
dynamics of any given time series is rarely of interest. An alternate approach, 
formulated entirely in terms of the statistical quantities of interest, is to attempt 
to solve the equations and boundary conditions which determine the Hopf char- 
acteristic function Z(J) (see [TJ [2]. Recently this approach has been applied to 
problems in atmospheric dynamics and stellar physics [3J d] ) . 

For a dynamical system with phase space variables X\ , • • • , X n , the characteris- 
tic function Z{J) is defined [H[2] as the time average of exp (i J ■ X(t)\ . Assuming 
ergodic behavior, 

Z(J) = ^exp (if ■ A(i)) } = J dn(X) exp (if ■ x) , (1.1) 

where dp(X) is the formal representation of the invariant measure. For a fractal 
which is not derived from a dynamical system, the characteristic function is defined 
solely in terms of a choice of measure dp(X). The small J behavior of Z carries 
information about the moments (Aj • • • Xj), via the Taylor-Maclaurin expansion. 
Remarkably, and less well known, is that the large J behavior carries information 
about the dimensionality of the chaotic orbit or fractal. We will be concerned with 
the large J behavior in this article. 

When the invariant measure can be written in terms of a probability density 
function, dp(X) = p(X)d n X, where p(x) is everywhere finite, Z(J) and p(X) 
are Fourier transforms of each other. However in many systems of interest, the 
dynamics is dissipative and a probability density function which is everywhere 
finite does not exist, although it may exist as a distribution. If an everywhere finite 
p did exist, it would satisfy -£ = —pV ■ v > along any trajectory, where v(X) is 

the velocity which is inconsistent with Poincare recurrence. In the presence 
of dissipation, the dimensionality d of the chaotic orbit is less than dimension n of 
the space spanned by A, and may in fact be non-integer. For a dissipative system 
with d < n, the Fourier transform of Z( J) can not converge. A sufficient condition 
for the existence of a Fourier transform is that Z( J) be C 2 , such that the integral 
J d n J Z (J)* Z (J) converges. Thus it would seem natural to define a dimensionality 
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corresponding to the maximum value of s < n such that the integral 

d n J\J\ s - n \Z(J)\ 2 (1.2) 

J|>£ 

converges. 

The characteristic function can be given exactly for some systems with integer 
dimension d, such as the Orszag-Mclaughlin model [11 j . In this and any other 
case in which the invariant measure can be written as p{X)d n X where p is product 
of finite functions with delta functions, it is not hard to see that the maximum 
value of s of which I s converges is indeed equal to the dimension of the chaotic 
trajectory; s max = d. The asymptotics of the characteristic function is not so 
trivially determined when the dimension is fractional. In fact, Z(J) need not fall 
off at large J, except in some averaged sense [H [7]. It is known that the maximum 
value of s for which I s converges is bounded above by the Haussdorf dimension d. 
I s is proportional to the 's-dimensional energy of the measure', defined by 

w = jj44W) i (L3) 

the inverse of which is known as the s-capacity. For all measures on a Borel set 
of dimension d in a complete metric space, the s-capacity vanishes for s > d. For 
s < d there exists a Borel measure such that the s-capacity is finite [8J- In this 
article we explicitly evaluate the asymptotic behavior of the Hopf characteristic 
function, or convergence criteria for I s , for some well known chaotic dynamical 
systems and fractals. 

The relation between fractal measures and the asymptotics of their Fourier- 
Stieltjes transform has also been discussed in [9j, where the asymptotics was char- 
acterized by a similar, although not manifestly equivalent, parameter. The context 
of that work was different, involving relating fractal spectral measures in quan- 
tum mechanics to the large time asymptotics of of temporal correlation functions. 
Translating to the present context, it was shown that 

- f J dJ'\Z{J')\ 2 ~ J~ v ' 2 (1.4) 
J Jo 

in the limit of large J, where T>2 is the correlation dimension of the fractal mea- 
sure. For the parameter s max which we have chosen to characterize the asymptotic 
behavior, the correspondence with another definition of dimension is unknown 
although, as noted above, s max is bounded by the Haussdorf dimension. 

For a simple choice of measures on the Cantor set, and generalized Cantor set, 
we will show that the bound s max < d is exactly saturated, so that s max = d. 
The convergence of I s for s < d is realized in non-trivial way, since Z{J) does 
not fall off with large J. Indeed, there is an infinite and unbounded sequence of 
values of J having the same positive Z(J). Numerical results for the characteristic 
function of the Lorenz attractor yields an estimate for s max , which is consistent 
with the known Hausdorff dimension. It must be emphasized that this estimate 
is far too crude to claim equivalence with the Hausdorff dimension, as opposed to 
the slightly lower correlation dimension. 
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2 Hopf characteristic functions for integer di- 
mension 



At present, chaotic dynamical systems for which the characteristic function is 
known exactly have integer dimension d, and the invariant measure can be written 
as products of finite functions with delta functions. In this case it is not difficult 
to see that the maximum value of s such that (|1.2|) converges is s max = d. 

An example of an exactly calculable characteristic function is that associated 
with the Orszag-Mclaughlin dynamical system |10| . defined by the equations of 
motion, 

dx ' 

— = X i+1 X i+2 + Xi-lXi-2 ~ 2x i+1 Xi-l (2.5) 

with i = 1, • • • ,n and periodic identification Xi +n = Xi. The Hopf characteristic 
function, computed in [TTJ, is given by, 

Z(J) = T{n/2){R\ Jl/2) 1 -™/ 2 J n/2 _M J\) , (2.6) 

where R is a constant related to a conserved quantity of the motion, and J n /2-i 
is a Bessel function. This characteristic function corresponds to a probability 
distribution 

P(x) = ^—5(\x\-R), (2.7) 
so that the dimension of the chaotic orbit is d = n — 1. For large | J|, 

l-n/2 



, , . / R\J\ \ 2 / , (n/2 - l)vr vr 

Z(J)«r(n/2) Ui J-^cos i?|J' 



(2.8) 

The large \ J\ scaling, \Z\ 2 ~ | J| 1_ ", implies convergence of (|1.2p for s max = n — 1, 
in agreement with the dimension d = n — 1. 



3 Hopf characteristic function of the middle 
third Cantor set 

The middle third Cantor set is a fractal defined by starting with the unit interval 
and removing the middle third, and then the middle third of each remaining seg- 
ment, ad infinitum. The box counting dimension D is defined by taking segments 
of length epsilon, and determining how the number of such segments necessary to 
cover the Cantor set scales as e — > 0. Taking e = (l/3) n , it is not hard to see that 
the number N required to cover the set is 2 n . N scales like (l/e) D with 

ln(3) 
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In this case the box counting dimension is equal to other common definitions of 
dimension, such as the Hausdorff dimension. 

To compute the characteristic function (exp(iJX)), we define a measure on 
the Cantor set by the n —> oo limit of a uniform probability distribution, with 
constant probability per unit length, on the set obtained after middle thirds have 
been removed n times. Writing the measure at each step as dp(x) = p n (x)dx, one 
has 

Pn(x) = (^j\ (3.9) 

such that J s dxp n (x) = 1. Here S n is the set obtained by removing the middle 
thirds n times. The characteristic function is 



lim ZJJ) (3.10) 

n— >oo 

where 

Z n (J) = / dxp n (x) exp(i Jx) (3-11) 

While the n — > oo limit of p n does not exist, the n — >• oo limit of Z n converges. 
Indeed, rapid convergence to a continuous function of is observed numerically. 
From the construction of the Cantor set, it is not hard to demonstrate that 

Z n+1 (3J) = ^(l + e 2iJ )Z n (J). (3.12) 

Existence of the limit (|3.1U|) and the scaling relation (|3.12|) implies 

Z(3J) = X - (1 + e 2iJ ) Z{ J) (3.13) 

Note Z(3 m irn) is independent of to, where m and n are is integers. One can check 
numerically that J = nn is not necessarily a zero of Z . Thus, it is not true that 
limj^oo Z{J) = 0, which would hold if the invariant measure could be written in 
terms of a finite probability density function. Perhaps contrary to expectation, 
Z(J) does not obey any simple large J asymptotics, such as a power law fall off. 
The scaling relation f)3. 12|) implies 

9J i-3J 

dj\Z(j)\ 2 r a = 3 1 - a / dj\Z(3j)\ 2 r a 

3J J J 

i r 3j 

= 2 3 Jj #(l + cos(2i))TO)rr a , (3 - 14) 

for any J. Thus, 



m\ 2 r a = Ysj 3 d i [—) xn(j)\zm 2 r a , (3.i5) 



>3J 
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where 



Xn = J] (l + coslS^^j)) • (3.16) 

m=l 

Note that Xn(j) does not converge to any function in the large n limit. In par- 
ticular, it diverges for any j = |£ with integer I and positive integer k. These 
values of j are a dense subset of the reals. Figure [T] shows X8,X9 an d Xw i n the 
neighborhood of j = ir/9. Despite the lack of convergence to any function, one 
can show that the n — > oo limit is a distribution, such that the limit 



lim / dj Xn (j)f(j) = (xJ) (3.17) 

exists for any function f(j) with a uniformly convergent Fourier series expansion^ 
on the interval [J, 3J]. The proof is given in appendix A. Rapid numerical con- 
vergence of (Xn, f) for some arbitrarily chosen test functions is shown in table 1. 
Assuming that lim n _ i>00 Xn(j) converges to a distribution sufficiently fast, then the 
sum (|3.15p behaves as 

n ^ ' 

for large n. Convergence of this sum requires 3 2 < 1. Therefore the maximum 
value of s such that 111. 2ft converges^], is 

ln(2) 

W = 1 - a = -)-(, 3.19 
ln(3) 

in agreement with the Hausdorff dimension. 



4 Hopf characteristic function of the gener- 
alized Cantor set 

A correspondence between the asymptotic behavior of the characteristic function 
and the fractional dimension also holds for the generalized Cantor set. The gen- 
eralized Cantor set can be defined as follows. Consider the map 

x -^f2r\x for x < — 

x ->1 - 2r](l - x) for x > - (4.20) 

1 We do not have a proof of convergence when the Fourier series is not uniformly con- 
vergent. 

2 One could also have attempted to apply the scaling relation Ij3.12j) to determine con- 

3 K j 

vergence of J Q dJ \Z{J)\ 2 J~ a for integer K as K — > oo. However this approach re- 
quires considerable care, as it involves evaluating {xk, f) with an integration over a region 
bounded by J = 0, where / = \Z{J\ 2 J~ a has a singularity. 
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n 


(Xn, fl) 


(Xn, h) 


1 


0.36677117 


-0.05595122 


2 


0.43466270 


-0.12962766 


3 


0.43529857 


-0.02549191 


4 


0.43545871 


-0.03928547 


5 


0.43546377 


-0.03951139 


6 


0.43546202 


-0.03948616 


7 


0.43546384 


-0.03948616 


8 


0.43546442 


-0.03986482 


9 


0.43546424 


-0.03987476 


10 


0.43546423 


-0.03987374 


11 


0.43546423 


-0.03987312 


12 


0.43546423 


-0.03987311 


13 


0.43546423 


-0.03987305 


14 


0.43546423 


-0.03987306 



Table 1: The integral f ± dJ Xn(J)f(J) evaluated for arbitrarily chosen test 
functions /i(J) = exp (— 5( J — 2) 2 ) and / 2 = cos(20J). Note the rapid 
convergence with increasing n. 



450 r 

400 

350 




0.3488 0.3489 0.3489 0.349 0.349 0.3491 0.3492 0.3492 0.3493 0.3493 



Figure 1: Plots of x&U) (blue), Xsij) (red) and Xio(j) (green) in the neigh- 
borhood of j = 7r/9 
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n 


(X?7,n> /l) 


{Xv,n, fy 


1 


0.07620179 


-0.04757256 


2 


0.11068899 


-0.01065505 


3 


0.08751986 


-0.00348152 


4 


0.08803434 


-0.00178861 


5 


0.08802082 


-0.00481686 


6 


0.08802103 


-0.00479188 


7 


0.08802102 


-0.00479230 


8 


0.08802103 


-0.00479224 


9 


0.08802103 


-0.00479226 


10 


0.08802103 


-0.00479226 


11 


0.08802103 


-0.00479226 


12 


0.08802103 


-0.00479226 


13 


0.08802103 


-0.00479226 


14 


0.08802103 


-0.00479226 



Table 2: The integral L v dJ Xv,n(J) J '(J) evaluated for r] = 5/4 and test 
functions /i(J) = exp (— 5( J — 2) 2 ) and f 2 = cos(20J). Note the rapid 
convergence with increasing n. 



for rj > 1. With each application of the map, points between ^- and 1 — ^- are 
mapped outside the unit interval [0, 1] while the segments [0, and [1 — ^,1] are 
stretched to fill the interval [0.1]. The generalized Cantor set is the set of initial 
conditions on the interval x = [0, 1] which are not mapped outside the interval 
after any number of applications of (|4.2(jp . For rj = | , this set is the middle third 
Cantor set discussed above. 

The characteristic function of the generalized Cantor set can be studied in 
essentially the same manner as the middle third Cantor set. Starting with the 
distribution p = 1 on the unit interval, a new distribution can be defined with each 
application of the map (j4.20p . With each application of the map, the set of initial 
conditions which remain within the unit interval is reduced, and the distribution is 
taken to be uniform (constant) on the remaining set of inital conditions, integrating 
to 1. For instance, after the first application of the map (|4.20p . the surviving initial 
conditions are the segments [0, ^] and [1 — -^, 1], on which the distribution is p = r/. 
After n applications of the map, one has 

Pn = if (4-21) 

on 2 n segments S n within the unit interval. The characteristic function is 

Z n (J)= / dxp n exp(i Jx) (4.22) 

JSn 

The definition of S n and p n yields the relation, 

Z n+1 (2 V J) = \ (l + e l2r>J{l -^ ] ) Z n (J) . (4.23) 
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Therefore, existence of the limit Z{J) = lim n _ 5 . 00 Z n (J) implies; 

Z(2i 1 J) = i (l + e i2riJ{1 ~^ Z(J) (4.24) 

Thus, 

dj \Z(j)\ 2 r a = (2r ? ) 1 " a / dj \Z(2 V j)\ 2 r a 

2r?J J J 

i r^ J / / 1 \\ 



2 (2vY~ a ]_ dj (1 + 008(2^(1--)) /)r.r". (1.2o) 



for any J. Therefore 

r°° r 2r i' J /f?n]l-«\" 

dj \Z{j)\ 2 r a = *" ) l^(i)| 2 J~ a , (4-26) 



where 



^[ (l + cos ((27?)— ^7(1 " ^))) • (4-27) 



The n — > oo limit of Xnn does not converge to a continuous function. However, 
we conjecture that it converges sufficiently rapidly to a distribution, such that the 
sum (|4.26p converges when — < 1. For several values of rj and for several 
arbitrarily chosen test functions /, we have numerically checked that (Xr],n,f) 
converges rapidly with increasing n, as illustrated in table 2. The maximum value 
of s such that (|1.2p converges, is then 

*- = 1 - a= ^' (4 - 28) 

which is equal to the Hausdorff dimension. 



5 Hopf characteristic function of the Lorenz 
attractor 

The Cantor sets considered above are self similar fractal sets, for which scaling 
arguments yield information about the asymptotic behavior of the characteristic 
function. The chaotic invariant sets associated with continuous dynamical systems 
are not generally self similar, making it difficult to determine the large J behavior 
of Z(J) analytically. In the following we give crude estimates for the asymptotic 
behavior of the Lorenz attractor, i.e. for the maximum value of s such that (jl.2p 
converges, by a numerical computation of Z{J). This approach is hampered by 
the fact that one can only compute Z up to a finite value of J, limited in size by 
the available computational power. 
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The Lorenz attractor [12] is defined by 

dx dy dz 

— = a(y -x), — = x(p-z)-y, — = xy - (3z. (5.29) 

We choose the value a = 10, p = 28, (3 = 8/3, and a set of N distinct initial 
conditions within the basin of attraction. These are evolved for a fixed time T. 
The characteristic function is given by 

TV t 

Z(J)= lim = -V- / dte if - Si ^ (5.30) 
i=l 







The choice of N is irrelevant; it may just as well be set to one. However it is 
computationally most efficient to choose a large N and evolve each set of initial 
conditions in parallel for a finite time to estimate Z{J). As the value of |J| is 
increased, one must increase N and/or the density of sampling of x(t) in time to 
avoid a problem akin to aliasing^. 

We wish to determine the minimum value of a = 3 — s such that the integral, 

1 = J d 3 J\J\- a \Z(J)\ 2 (5.31) 

converges Because of computational constraints, we will instead pick a three- 
dimensional unit vector n and demand convergence of the integral 

/>oo 

I{n)= dJ J 2 J- a \Z(Jh)\ 2 (5.32) 
Jo 

For the few choices of h we have made, the result seems to be insensitive to the 
choice of h. A contour plot of the integral (|5.32j) as a function of a and the endpoint 
of integration J en( j is shown in figure [21 suggesting a critical a in the neighborhood 
of a = 1, above which the integral converges and below which it diverges. This 
is not inconsistent with saturation of the bound on the dimension, d > 3 — a cr it, 
since the Hausdorff dimension (see [15] for a recent estimate) is d ~ 2.06. 

Of course, numerical results are extremely crude and merely suggestive of 
asymptotic properties, since the they are limited to finite J endpoints of integra- 
tion. A more clever approach is required to accurately determine the asymptotics 
of Z(J). It would be very interesting if there existed some asymptotic analysis of 
the differential equations which the characteristic function satisfies (see e.g. |llj). 
which could determine a cr it- As in the case of the Cantor set, there is no reason to 
expect the large J behavior of Z{J) to a have a simple functional form, such as a 
power law fall off. Thus it may be that one must consider some suitably averaged, 
or integrated form of these differential equations to obtain a C rit- 



6 Conclusions 

While the small J behavior of the Hopf characteristic function Z(J) is related to 
correlation functions, via the coefficients of the Taylor-Mclaurin expansion, the 



3 The rate at which sampling densities must be adjusted with increasing | J| is presum- 
ably related to the box counting dimension. 
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a 




100 200 300 400 500 600 700 800 900 

end 

Figure 2: Contour plot of the integral jf cnd dj j 2 ~ a \Z(jz)\ 2 as a function of 
a and J cn d- 

large J behavior is related to dimensionality. For the parameter we have chosen to 
measure asymptotics, the large J behavior of the characteristic function satisfies a 
bound given by the Hausdorff dimension of the set. We showed this bound to be 
saturated for a natural measure on the generalized Cantor set. Crude numerical 
results are also consistent with saturation for the Lorenz attractor. A number of 
interesting questions remain un-answered. 

In order to make the arguments relating the dimensionality of the Cantor set 
to the asymptotic behavior of Z{J) more rigorous, it remains only to prove that 
that the functions Xn, defined in (|3.16p and (|4.27[) approach a limiting distribu- 
tion sufficiently fast at large n. Thus far we have only shown (see appendix A) 
that linin^oo Xn is a distribution on the space of functions which have absolutely 
convergent Fourier series expansions. A better understanding of the properties of 
the functions Xn would be valuable. Note that they do bear some resemblance to 
Weierstrass's non-differentiable function. 

Subject to certain assumptions about analyticity, the small J behavior of Z{J) 
is related to the large J behavior. In principle, the large J behavior could be 
determined by a re-summation of the cumulant expansion ln(Z) = Y^=o c nJ n ■ 
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The dimensionality should give constraints on the behavior of the cumulants c n 
for large n. It could be very interesting to elucidate these constraints. 

Perhaps the differential equations satisfied by the Hopf characteristic function 
of a chaotic dynamical system, which follow from the equations of motion, can 
somehow be used to draw conclusions about the large J asymptotics. If so, such 
an approach will likely involve an integrated or suitably averaged form of these 
equations, since Z(J) need not display a simple asymptotic fall off when the di- 
mension is non-integer. Indeed, we have seen that there is no asymptotic fall off 
for the Cantor set, nor do we observe a power law fall off for the Lorenz attractor. 

This work was motivated in part by [13} 114] . in which a class of chaotic dy- 
namical systems was found for which exact statistics, in the form of a probability 
distribution function, is known. These dynamical systems were reverse engineered 
starting with a probability distribution and a two-form. Due to the existence of 
a finite probability distribution function, these systems do not have fractional di- 
mension, and are dissipative in some regions of the chaotic orbit but not globally. 
There may be a way to extend the inverse approach to dissipative systems with 
fractional dimensions, starting with a characteristic function having no Fourier 
transform, rather than starting from a probability distribution. 
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Appendix A 

Here we show that the n — > oo limit of 

n 

Xn= [J (l + cos^-^j)) . (A.33) 

m=l 

is a distribution, such that the limit 

lim / djxn(j)fU) (A.34) 

exists for any function f(j) with a uniformly convergent Fourier series expansion 
on the interval [J, 3 J]. 
In light of the identity 

cos(a + b) + cos(a - b) 
cos(a) cos(o) = , (A. 35) 

the highest order term in the Fourier expansion of \n is cos(pj) with 

n 

p = 2 ^ = 3"- 1 (A.36) 

m=l 
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It follows from (|A,35h that the Fourier expansion of Xn+i = (1 + cos(3 n 2j))xn will 
have the same coefficients as that of Xn for all modes of order p or less, with higher 
order modes having coefficients less than or equal to the maximum coefficient of 
Xn, which by induction is 1. Therefore Xn has a Fourier expansion in cosines of 
even integer multiples of j with positive coefficients less than or equal to one. 

For simplicity, let us now choose J = ir/2 and exchange f(j) in (|A.34j) with 
a 7r periodic function which is equivalent to f(j) in the integration range. Now 
f p (j) and Xn have Fourier expansions 



f P d) = Yl cos ( 2k i) + h sin(2/cj)) , 

k=0 
3 n -l 

Xn(j) =Y, Ck cos ( 2k 3)> < c fe < 1, (A.37) 



k=0 

such that 

"371-/2 

If the Fourier series expansion of fp(j) is absolutely convergent, i.e. 



/ dj f P (j)Xn(j) = vra c + 77^ afcCfe - (A.38) 

Jn/2 1 ,. -, 



^2\a k \<oo (A. 39) 



k=0 

converges, then by the comparison test 

3 n -l 



lim V \a k c k \ (A.40) 



n— >oo 

k=0 



must also converge. Thus (|A.38|) converges in the limit n — > 00, and we conclude 
that lin^^oo Xn is a distribution. 
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